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Based on the combined model of the double exchange and the polaron, we have studied the small-to- 
large polaron crossover transition and explored its effects on the magnetic and transport properties 
in colossal magnetoresistance (CMR) manganites. We have used the variational Lang-Firsov canon- 
ical transformation, and shown that the magnetic and transport properties of both high and low Tc 
manganites are well described in terms of a single formalism. We have reproduced the rapid resistiv- 
ity drop below Tc, a realistic CMR ratio, and the first-order-like sharp magnetic phase transition, 
which are observed in low Tc manganites. 
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Despite intensive researches on the colossal magnetore- 
sistance (CMR) manganites Ri-^Aj^MnOa (R = rare- 
earth; A = divalent cation), there are many unusual prop- 
erties that remain to be clarified. The most outstand- 
ing feature of these compounds is the close relationship 
between transport phenomena and the magnetism for 
0.2 ^ X ^ 0.5, which have been traditionally understood 
by the double exchange (DE) mechanism ||l|,||. However, 
many recent experiments, especially the giant isotope ef- 
fect in the magnetic transition temperature (Tc) B, as 
well as theoretical works H] indicate that the electron- 
phonon interaction plays a crucial role in the physics of 
the CMR manganites. The small-to-large polaron transi- 
tion near Tc has been reported by various experimental 
tools tH. 

The Ro.yAo.sMnOa CMR compounds have a wide 
range of Tc from 70 K to 365 K depending on the av- 
erage radius of R and A cations {{ta)) 0- Interest- 
ingly, the low Tc manganites are different from the high 
Tc manganites in many properties. The low Tc man- 
ganites have higher resistivity, much sharper resistivity 
peaks near Tc, and more enhanced magnetoresistance 
(MR). They often exhibit hysteresis behaviors in the 
temperature dependent resistivity [^. Further, in con- 
trast to the other manganites that experience a metal- 
insulator (M-I) transition near Tc, the highest Tc man- 
ganite Lao.rSro.aMnOa is metallic in the whole tempera- 
ture range. 

The magnetic mo- 

ment measurement for Lao.67(Baa;Cai_a;)o.33Mn03 with 
varying x shows that the magnetic transition becomes 
steeper for low Tc samples, reminiscent of the first-order- 
like phase transition for x = 0.0 (Tc=340 K and 265 K 
for X = 1.0 and x = 0.0, respectively) [Q. Also the 
NMR measurement for Pro.TCao.isSro.igMnOg (rc=180 
K) and Pro.rBao.sMnOa ((Tc=175 K) shows that the hy- 
perfine fields on ^^Mn nuclei remain finite up to Tc, in- 
dicating the first order magnetic transition [|2[ . Neutron 
scattering experiment indicates that the spin dynamics 



near Tc of low Tc Ndo.ySro.aMnOa (Tc=197.9 K) is very 
distinct from that of high Tc Pro.rSro.sMnOs (Tc=300.9 
K) ||l^ . The spin wave stiffness data for the former show 
no evidence of the expected spin wave collapse at Tc, 
as in Lao.eyCao.saMnOa |1J]. It is pointed out that the 
exotic spin dynamics in the low Tc manganite is related 
to the increased electron- lattice coupling |13] . 

Various theoretical models are proposed to explore the 
origin of the CMR phenomena. Among those, the com- 
bined model of the DE and the polaron explains the 
transport and other many experimental features rather 
well |l5|-[l7||. Indeed, this model could explain the resis- 
tivity and the magnetic behaviors of the high Tc man- 
ganites such as Lao.TBao.aMnOa. However, it appears 
that this model cannot describe a rapid drop of the re- 
sistivity below Tc and very large MR ratio of the low 
Tc manganites such as Lao.rCao.aMnOa. Including the 
effect of the phonon hardening below Tc improves the 
resistivity drop and the MR ratio substantially, but the 
characteristic magnetic field needed to reproduce the ob- 
served MR is still too high {^ 15 — 20 T), as compared 
to the experimental value (^ 4 T) |1^. It is thus ar- 
gued that the polaron model is not sufficient to explain 
the physics of CMR |jl9]. Therefore the central issue is 
whether the other mechanism than the the DE and the 
polaron should be invoked to understand the physics of 
CMR manganites. 

Motivated by the very distinct behaviors in the resis- 
tivity and the magnetism between the low and high Tc 
manganites, we have reexamined the effects of the po- 
laron transition in the combined model of the DE and the 
polaron. Employing the variational Lang-Firsov (VLF) 
canonical transformation, we have investigated the small- 
to-large polaron transition and explored its effects on 
magnetic and transport properties of both the low and 
high Tc manganites. Roder et al. ||l^ have tried a similar 
approach, but they did not show the transport proper- 
ties explicitly in connection with the CMR phenomena. 
Millis et al. Ilq] have obtained comparable results using 



the dynamical mean field approximation, but the physics 
producing CMR is hidden and the resulting characteris- 
tic magnetic field needed to reproduce the observed MR 
is also too high. 

In order to study the polaron transport incorporating 
the DE, we adopt the following Hamiltonian of the Zcncr- 
Holstein type: 



H = -t-i{T) Y^ 4^^-c, + Y^ ujqa\aq 



iS 



E4^ 
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Mg{aq + alg). 
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Here the bandwidth variation by the DE is included in 
the hopping parameter with a factor 7(T) — (cos(6'/2)), 
where 9 is the angle between local spins of neighboring 
Mn sites. Mq is the parameter representing the strength 
of the electron-phonon interaction. For simplicity but 
without loss of the physics, we consider the single e^ or- 
bitals and the single phonon modes, which is plausible to 
describe the metallic phase of CMR manganitcs [|l^ . 

In the strong coupling {Ep = ^ M"^ /uOq ^ t) and the 
non-adiabatic limit [uq ^ t), one can apply the conven- 
tional Lang-Firsov (LF) transformation and obtain the 
resistivity via the Kubo formula, considering the hop- 
ping term as a perturbation [ [L7[ . In the opposite weak 
coupling limit, the analytic solution can be obtained by 
treating the interaction term as a perturbation. It is, 
however, difficult to evaluate physical quantities in the 
intermediate regime inbetween. 

We have adopted the VLF transformation to describe 
the small-to-large polaron transition and the associated 
physics in the intermediate regime. In the VLF method, 
one performs the Lang-Firsov transformation with some 
variational parameters and finds their values from the 
minimization conditions of the total energy. We follow 
the method of Das and Sil (20| in which three parameters 
are assumed: on-site lattice distortion A, nearest neigh- 
bor distortion A', and two-phonon coherent state param- 
eter a. The coherent state allows for the anharmonicity 
of the lattice fluctuations, which become important at 
finite polaron densities and in the intermediate coupling 
and frequency regime pl[ |. The nearest neighbor dis- 
tortion can be neglected in the three dimensional (3D) 
case for which the small-to-large polaron transition is so 
abrupt that A' does not contribute much in the course 
of the transition. 

Using the VLF transformation, we have 
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phonon coherent state |'0ph) — e "^q^"'''' "> ""j" i' 
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(a > 0) as a trial phonon wave function, one can obtain 
the effective polaron Hamiltonian. Now the ground state 
energy is calculated in the framework of the mean-field 
theory: 



E/N = -2zt-f{T) exp 
~nEp{2 ~ A)A 



.Ely 



n(l — n) + ujQ sinh (2a) 
(3) 



where F = A^e~^", z is the coordination number, and 
n — N^^^^rii is the carrier density. We assumed a 
single phonon frequency wq. The values of A and a are 
obtained from the minimization of Eq. (H). With A = 1.0 
and a = 0, that is, with F = 1.0, one can reproduce the 
result of the conventional LF transformation. Note that, 
in the VLF formalism, the effective hopping parameter 

is given by t = i7(r) expf — -^F j with an additional 

parameter F(< 1.0) in the polaron narrowing factor, so 
that the narrowing effect becomes weakened in the VLF 
formalism. 

Figure provides F and the polaron narrowing fac- 
tor as a function of electron-phonon coupling constant 
Xp = Ep/{2zt) for both the adiabatic and non-adiabatic 
cases. As shown in Fig. ||, the large-to-small polaron 
transition, the so-called self-trapping transition takes 
place near Ap « 1 and is steeper for the adiabatic case. 
Also notable is that the critical value of AJJ at which 
the transition occurs is smaller for the adiabatic case. 
Though we do not show the results for ID and 2D sys- 
tems, the transition for 3D systems is steeper. These 
features are consistent with those obtained using other 
methods such as the exact diagonalization [p2|, the dy- 
namical mean field approximation p3[ , and the quantum 
Monte-Carlo method [|§|2|]. Since the CMR manganites 
are close to 3D adiabatic systems, a very steep small-to- 
large polaron transition is expected. 

It is tempting to interpret that the very steep self- 
trapping transition in Fig. |l| is the real first order phase 
transition. Caution is needed, however. According to the 
study of Lowen [g6[ for the Hamiltonian of the Holstcin 
type, the self-trapping transition should be an analytical 
crossover not an abrupt (nonanalytical) phase transition. 
Indeed, numerical calculations for various model systems 
p2[- p5tp7[ show rather continuous self-trapping transi- 
tions. Then the seemingly nonanalytical transition seen 
in Fig. n may not be real but result from the shortcom- 
ings of the VLF method itself. Nevertheless, the general 
trend of the self-trapping transition obtained in Fig. m can 
be applied to CMR manganites for the analysis of trans- 
port and magnetic properties, by using a well-behaved 
function that smooth the nonanalytical transition (dots 
in the bottom panel of Fig. p^(a)). 

To extend Eq. (|3|) to the finite temperature regime, 
the entropy term originating from the magnetic ordering 
must be included |E8[. The resulting free energy per site 
is represented by 
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-2ztj{i') exp 
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A)A - ksT {HQsiiy)] - vms{v)] , (4) 
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rameter, Qs{iy) = X]m=-s ^'^pC^-^/'^) is the partition 
function with spin 5 = 2, and ms{v) is the normahzed 
magnetization per site. A^o is the phonon distribution 
function and can be ignored in the temperature region of 
interest. ^{T) is also a function of ^', being constant (0.6) 
for T > Tc and increasing below Tc up to 0.8. Now the 
free energy is to be minimized instead of Eq. (ra) . 

Upon cooling below Tc, the DE hopping parameter 
tj{T) increases. Accordingly, if Xp decreases below Ap, T 
becomes reduced a lot to weaken the polaron narrowing 
effect as seen in Fig. ^. This implies that F is an implicit 
function of the hopping parameter: with a larger hopping 
having reduced F. In consequence, the effective hopping 
parameter t is rapidly enhanced to cause an enormous 
reduction of the resistivity below Tc (see Fig pi) . In turn, 
this has an effect of increasing 7(T) to make the magnetic 
transition much steeper than in the case of the DE onl 
model. This feature will be discussed more in Fig 
The external magnetic field gives rise to the same effect, 
producing the CMR along with the bandwidth increment 
(Fig. I). 

From Eq. (H), Tc ior S — 2 can be derived: 
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This relation provides a possible account of the intriguing 
evolution of Tc with chemical pressure in Ro.rAo.sMnOs, 
determined by (r^) flQ]. This behavior is qualitatively 
understood in terms of the variation of the {Mn — O) 
bond length with (r^) and the subsequent variation of 
the hopping parameter t p9| |. However, the observed 
variation in the bond length, at most 1%, is thought to 
be too small to explain ~ 500% variation in Tc- However, 
once employing the above polaron narrowing effects, one 
can account for the huge variation of Tc even with the 
small variation of the hopping parameter t. That is, the 
increase of (r^) enhances i, and so reduces Ap. As a 
result, F in the exponent of the polaron narrowing factor 
in Eq. (||) decreases so as to enhance Tc a lot. In a 
similar way, the decreasing isotope effect with (r^) can 
also be understood by noting that the isotope exponent 
/3(= — AlnTc/AlnMo, Mq'- oxygen isotope mass) is a 
decreasing function of the hopping parameter through F: 

/3~^r(i)§. 

Since we have determined the polaron narrowing fac- 
tor for all the range of Xn (Fig. |l|), Apo for given Tc can 
be estimated using Eq. (g) . We have larger Apo for lower 
Tc to have a higher slope of F at Apo, e.g., Apo = 0.646 
and 0.619 for Tc = 100 K and 400 K, respectively, for 
parameters given in Fig. 0. Hence, as 7(T) increases 
below Tc, the subsequent variation of Ap drives a large 



reduction of F. Therefore one can expect a steeper mag- 
netic transition and more rapid drop of the resistivity 
for low Tc manganites. Figure presents the tempera- 
ture dependent behaviors of the magnetization for vari- 
ous Tc- It is seen that the magnetic transition becomes 
steeper for lower Tc- Thus, the magnetic transition for 
Tc < 200 K is almost like the first order phase tran- 
sition. Note, by contrast, that the magnetic transition 
in the DE only model takes place much more smoothly 
refiecting the second order phase transition. In view 
of these features, the change of magnetic properties in 
Lao.67(Ba:rCai_a;)o.33Mn03 with varying x can be eluci- 
dated ||ll[]. Further, this feature provides a clue to de- 
scribe the exotic spin dynamics observed in low Tc man- 
ganites. Since the transition is close to the first order 
phase transition, the coexistence between the paramag- 
netic and the ferromagnetic phase can occur in the course 
of the transition, manifesting features of the discontinu- 
ous drop of the hyperfine fields |12| and the spin wave 
stiffness |13] near Tc- 

The routine evaluating the conductivity by means of 
the Kubo formula is the same as in Ref. |Ol except for 
F in the polaron narrowing factor. FigureH shows the 
resulting resistivity as a function of temperature. Here 
the phonon hardening effect below Tc is also taken into 
account ||l^ . The M-I transition above Tc is clearly seen 
in the cases of low Tc- For Tc = 100 K, the resistiv- 
ity peak is very sharp and the MR is very large. Note 
that the y-axis is in the log scale so that the resistivity 
drop below Tc amounts to three orders of magnitude. 
Hence the characteristic field for the observed MR be- 
comes low enough, in agreement with the experiment. 
In the case of Tc = 200 K, we have ~ 94% MR ra- 
tio for H = 5 T. In addition, the lower Tc manganites 
show higher resistivity in consistent with the experiments 
|10| . For higher Tc, the peak becomes dull, and even- 
tually the M-I transition disappears for Tc larger than 
300 K. This is again in agreement with the observation 
for Lao.ySro.aMnOa. One must note that the magnetic 
transition (Fig. 0) and the resistivity drop (Fig. 0) are 
steeper for low Tc manganites, demonstrating that both 
the CMR transport phenomena and the magnetic transi- 
tion are closely correlated and are governed by the same 
physics. 

Finally, it should be pointed out that the important ef- 
fects, such as the double degeneracy of the Cg orbitals, the 
inter-orbital Coulomb correlation, and the Jahn- Teller ef- 
fects, are not considered in the present study. One can 
in principle incorporating these effects, following the for- 
malism by Zang et al [|30| . The Jahn- Teller interaction 
yields the effective DE besides the Zener DE, and the 
Coulomb correlation gives rise to an additional band nar- 
rowing effect. Therefore, these effects are expected to en- 
hance further the localization effect and the essential re- 
sults of the present study will not be altered. The study 
that takes into account these effects explicitly is under 
progress. 

In conclusion, we have investigated the small-to-large 



polaron crossover transition for the combined model of 
the polaron and the DE interaction and explored the ef- 
fects of the polaron narrowing on the magnetic and trans- 
port properties of both the low and high Tq CMR man- 
ganites. We have shown that the polaron narrowing ef- 
fect is much more pronounced in the low Tc manganites, 
which explains the rapid resistivity drop, large MR ratio, 
and the first- order-like sharp magnetic phase transition. 
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FIG. 1. The polaron narrowing factor exp 



and the 



n = 0.3 obtained using the VLF transformation: (a) the 
adiabatic case {ujo/t — 0.13) and (b) the non-adiabatic case 
{uJo/t = 3.0). Results of the conventional Lang-Firsov (LF) 
transformation (F — 1.0) are also presented in the top panel 
for comparison. Dots in the bottom panel of (a) represent an 
illustration of the continuous well-behaved function smooth- 
ing the abrupt self-trapping transition. 
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FIG. 3. Temperature dependent resistivity without and 
with the magnetic field of 5 T (lower curves) for various Tc's. 
Note that the j/-axis is in the log scale. Parameters used in 
the calculations are the same as in Fig. S. 
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FIG. 2. Magnetization vs. temperature for various Tc's. 
Result in the DE only model is also presented for comparison. 
Parameters used in the calculation are t'y{Tc) = 0.3 eV and 
ujo = 0.04 eV. 



